Double-layer shocks in a magnetized quantum plasma 
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The formation of small but finite amplitude electrostatic shocks in the propagation of quantum ion- 
acoustic waves (QIAWs) obliquely to an external magnetic field is reported in a quantum electron- 
positron-ion (e-p-i) plasma. Such shocks are seen to have double-layer (DL) structures composed 
of the compressive and accompanying rarefactive slow-wave fronts. Existence of such DL shocks 
depends critically on the quantum coupling parameter H associated with the Bohm potential and 
the positron to electron density ratio <5. The profiles may, however, steepen initially and reach a 
steady state with a number of solitary waves in front of the shocks. Such novel DL shocks could be a 
good candidate for particle acceleration in intense laser-solid density plasma interaction experiments 
as well as in compact astrophysical objects, e.g., magnetized white dwarfs. 
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Recent studies [IHl] have shown that the formation 
of stationary current free double layers (DLs) may be 
possible in dense plasma environments where the forces 
namely, (i) quantum statistical pressures for electrons 
and positrons, (ii) electron and positron tunneling asso- 
ciated with the Bohm potential play important roles in 
the propagation of quantum ion-acoustic waves (QIAWs). 
Inclusion of these quantum forces along with electron and 
positron angular momentum spin allows also the exis- 
tence of very high-frequency dispersive electrostatic and 
electromagnetic waves (e.g., in the hard X-ray and 7-ray 
regimes) with extremely short wavelengths (see for recent 
review in quantum plasmas, Ref. Q). However, there 
was no indication or explanation of the forward propa- 
gating (upstream) DLs that may exist, especially in a 
strongly magnetized quantum plasma, as a form of com- 
pressive and accompanying rarefactive slow-wave fronts. 
Formation of these upstream DLs is one of the most strik- 
ing features of the magnetohydrodynamicXMHD) shocks, 
and has been observed by Tajima et al [6j, possibly for 
the first time, in a classical MHD flow. Moreover, there 
is also a wide-spread interest in investigating DLs as a 
possible acceleration mechanism in various space and as- 
trophysical plasma environments (see, e.g., @-[Hj])- In 
a thought-provoking series of discussions, Alfven consid- 
ered DLs to be a central paradigm in plasma astrophysics 

Our current knowledge of DLs in plasma physics, how- 
ever, is insufficient for us to judge with much confidence 
and belief what roles DLs may, indeed, play especially in 
astrophysical environments. Thanks for the evidence of 
particle acceleration in a magnetized white dwarf which 
has been reviewed by Jager in the past [T3|. It has been 
reported there that if the current density is large enough 
in a tenuous magnetosphere, DLs can be formed leading 
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to a large electric field, and hence monoenergetic elec- 
trons and ions with higher energies. Recently, a new 
study (see, e.g., Ref. demonstrates that a white 

dwarf star may pulse like a pulsar, e.g., AE Aquarii can 
emit pulses of high energy X-rays as it rotates on its 
own axis. It has also been predicted that since pulsars 
are known to be sources of cosmic rays, white dwarfs 
should be quiet but numerous particle accelerators 
contributing many of the low-energy cosmic rays in our 
galaxy. One may thus be interested to know the physical 
mechanisms of such particle accelerations in these astro- 
physical compact objects. In this regard it may be noted 
that ever since the discovery of cosmic rays, the problem 
of understanding their origin and acceleration has not yet 
been fully understood. In this context, a number of mod- 
els and several processes have been proposed in the early 
days for their origin and the acceleration mechanism (see 
for some discussions, e.g., Ref. 16]). Furthermore, in 
addition to the shock-wave acceleration mechanism pro- 
posed early by Colgate and white 17|, Goldreich et al 
[l8j and Gunn et al [l9| had also shown that pulsars can 
accelerate particle to very high energies. Recently, the ex- 
pected abundance of cosmic-ray electrons and positrons 
from pulsars and magnetars have been studied [201 ] . 

A basic prerequisite for particle acceleration in pulsars 
is the presence of magnetospeheric plasma. On the basis 
of which Goldreich and Julien [l8| had shown that the 
large electric field due to a nonzero E.B at the neutron 
star surface, can overcome the gravitational potential and 
ensure a minimum plasma density. They also showed 
that the plasma flowing along the co-rotating magnetic 
field lines and escaping out at the light cylinder (defined 
by the distance from the center of the neutron star) will 
experience a potential drop and hence be accelerated. 
Furthermore, since an oblique magnetic rotator (in which 
the magnetic and rotational axes do not coincide) emits 
electromagnetic (EM) waves of the same frequency as 
that of rotation, particles can be accelerated at or beyond 
the light cylinder to high energies extremely efficiently by 
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the low-frequency EM waves emitted by a pulsar. This 
was first proposed by Gunn et al [l9| . 

Note, however, that in addition to degenerate elec- 
trons, there would also exist degenerate positrons, e.g., 
in magnetars and in the next generation of intense laser- 
solid density plasma interaction experiments [j| for which 
DLs could well be responsible for ion acceleration therein. 

Motivated by these facts, we report in this work, 
the existence of forward propagating slow-wave shocks 
in a magnetized quantum plasma, whose constituents 
are electrons, positrons and positive ions (hereafter re- 
ferred as e-p-i) . Because of the sufficient lifetime of the 
positrons compared to the ion time scale, plasma can 
become an admixture of electrons, positrons and ions. 
Such e-p-i plasmas are believed to exist, e.g., in the mag- 
netosphere of pulsars, in the active galactic nuclei, in 
the regions of the accretion disks surrounding the cen- 
tral black holes [H, [2ll-[23j. Van Allen radiation belts, 
near the polar cap of fast rotating neutron sta rs [24j . 
supernova remnants [25(, in intense laser fields 
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compact astrophysical objects (e.g., giant planetar 
teriors, white dwarfs, neutron stars/ magnetars) [2 
tokamaks [28] as well as in the early universe (29j . 
that the process of electron-positron (e-p) pair creation 
and annihilation may occur in relativistic plasmas at high 
temperatures, when the temperature of the plasma ex- 
ceeds the rest mass of electrons. However, for the prop- 
agation of ion-acoustic waves in e-p-i plasmas, the e-p 
pair annihilation can be neglected in the sense that the 
electron-positron lifetime is much larger than the charac- 
teristic time scale for collective oscillations (see for detail 
derivation and discussion, e.g., Ref. [30|). 

In what follows, we will consider the quantum force 
associated with the Bohm potential to provide higher or- 
der dispersion along with the charge separation effect as 
well as the magnetic-field-induced dispersion anisotropy. 
The ions are assumed to be cold and the motion is consid- 
ered on the ion-acoustic time scale. Because of their light 
masses, electrons and positrons will be highly magnetized 
compared to the ions ( ion Larmor radius is much larger 
than that of electrons or positrons), and will move almost 
parallel to the external magnetic field, so that electrons 
and positrons may be described by the quantum modified 
Boltzmann-like distributions. We will derive a modified 
Korteweg-de Vries (MKdV) equation with a quadratic as 
well as cubic nonlinearity that describes the dynamics of 
quantum ion-acoustic (QIA) DLs, and investigate some 
interesting properties of such DLs for large times. We 
observe a novel DL structure to the upstream shocks, 
which is composed of a compressive slow- wave shock and 
a rarefactive slow- wave front. 

Under the above assumptions, the motions of ions, 
electrons and positrons in the propagation of QIAWs 
obliquely to the external magnetic field B = Bqz can 
be described by the following equations |3l| . 

(1) 



dtrit + V.(njv) = 0, 
(d t + v.V)v= -V^ + w c i 



tt2 

n 2 J 3 = 1 + 20+ — = V 2 V^I, (3) 
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:V 2 ^, 



V 2 </> = n e — Sn p — (1 — 5) 



(4) 



(5) 



where n a is the number density of a— species particle nor- 
malized by the unperturbed value n Q o, v = (v x ,v y ,v z ) 
is the ion fluid velocity normalized by the ion-acoustic 
speed c s = ^Jk^Tj^Jmi with fc^ denoting the Boltz- 
mann constant, Tp e the electron Fermi temperature and 
rrii the ion mass. Also, </> is the electrostatic potential 



normalized by ksTp e /e, 



eB^/rriiUjpi is the ion- 



cyclotron frequency normalized by the ion plasma fre- 
quency, uipi = \JrtiQe 2 /eorrii, S — n p o / n e o is the positron 
to electron density ratio, H = hLu pe /kBTp e is the quan- 
tum parameter denoting the ratio of the electron plasmon 
energy density to the Fermi thermal energy. The space 
and time variables are normalized by the Fermi Debye 
length \p = c s /uj p i and the ion plasma period re- 
spectively. In Eqs. ([3|), (HI we have used the following 
Fermi-Diracpressure law due to electron and positron 
degeneracy [321 ] . 
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where Vp a = \Jk^Tp^Jrn a is the Fermi thermal speed. 
In Eqs. , , the terms proportional to H associated 
with the Bohm potential, account for typical quantum 
effects such as tunneling. In a broad sense, we refer to 
these particularities arising from the wave-like nature of 
the charge carriers as 'quantum diffraction effects'. 

In order to describe now the dynamics of propagat- 
ing DLs, we employ a reductive perturbation technique 
in which the independent variables are stretched as 
£ = e(lx% + lyU + lz% ~ Mt), t = e 3 t, where e is a 
small parameter representing the strength of the wave 
amplitude, M is the phase speed normalized by c s and 
l x ,l v ,l z are the direction cosines of the wave vector 
along the axes such that I 2 + I 2 + I 2 — 1. The de- 

oo 

pendent variables are expanded as n a = 1 + ^2 e J n a j, 

3=1 

OO oo 

Here the transverse velocity components (v x ,v y ) appear 
at higher order of e than that of the parallel component 
v z . This anisotropy is introduced by the influence of 
strong magnetic field. In this approximation the ion gy- 
romotion is treated as higher order effect. By inserting 
these expressions into Eqs. CO)-© and collecting the 
terms in different powers of e, we obtain in the lowest 
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where a = 8 1 / 3 . The flow may be outward or inward 
depending on the sign we consider in Eq. ([5]). Moreover, 
M < 1, i.e., the QIAWs propagate with the phase speed 
smaller than the ion-acoustic speed, and M increases as 
the positron to electron density ratio increases. 

Next, the coefficient of e 2 vanishes leading to Q4>\ = 
where Q = 1/a - 9(1 + a)/ {1-8). Without loss of 
generality we may assume that e<pi — > — </>o/2 as r = 
(x, y, z) — > oo so that 4>2 — > as r — > oo. Thus, Q should 



be at least of the order of e, i.e, 



e and is to 



be added to the third order contribution from Eq. (|5|). 
This gives a favorable condition for DL shocks instead of 
solitons. Physically, under this condition the free as well 
as the trapped particles are assumed to adjust themselves 
rapidly in order to maintain the quasineutrality at any 
time on each side of the propagating shocks. Now, the 
vanishing of the coefficients of e 3 gives five equations. 
When <9(/> 3 /<9£, dna/d£, dv z z/d£ etc. are eliminated and 
the first-order quantities are inserted into the resulting 
equation, terms containing 4>2 and </>3 cancel, and the 
following MKdV equation is obtained (0 = e(f>i). 



2 dcf) 



(8) 



where A = P/S, \i = Q/S,j = R/S with P = 6(l+a)/l z , 
R = 135(1 + a) 3 / (1 - S) 2 - 3 (1 - 8) /28 and 5 = 1 + 
(l-5)(l-Zf)/u^-9F 2 (l + a)/4a. An asymptotic shock 
solution of Eq. © can be obtained as [33[ 



>(£,t) = -ytanh 
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(9) 



where £ = £ — Vt, S\ is a constant and N — 1 represents 
the number of solitary waves in front of the shocks. Let 
us first investigate analytically the coefficients A, /i and 
7 with the system parameters. We see that for 5 < 1, 
P is always negative and Q(~ e) ^ according as < 
S < 0.001 or 0.001 < 8 < 1. Also, R is negative for 
< 8 < 0.0065 and positive otherwise. Moreover, 5^0 
according as H ^ H c , where H c is the critical value of H 
given by 
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(10) 



Note that this critical value, which depends paramet- 
rically on the density ratio 5, the obliqueness parame- 
ter l z and the ion-cyclotron frequency w c gives a criti- 
cal value of the electron density. The smaller the val- 
ues of H c the larger are the electron number densities. 



In order to consider the smaller values of H c or higher 
densities, one might have to disregard the charge separa- 
tion effect (the unity in the square brackets) in strongly 
magnetized plasmas (uj c < 1). In this case the criti- 
cal value scales as H c ~ (1 + o)/3. In a weakly 
magnetized case, H c ~ 2y / a(l — <5)(1 — Z 2 )/ (1 + a)/3oj e , 
which can relatively be larger than the strongly magne- 
tized case. However, in both the cases one has to be 
careful about the particle density range in which the 
Fermi thermal speed is much smaller than the speed of 
light in vacuum (Vp a << c) and the coupling parameter 
satisfies the relation: gQ = 2m e e 2 / £ah 2 (37r 2 y/no) 2 ^ 3 < 1 
(this corresponds to the case where the quantum collec- 
tive and mean field effects are important). Since Q ~ e 
for S ~ 0.001, we consider the regime < 6 < 0.0065 in 
which R < 0. Thus, in order that the DL solutions exist 
we must have H > H c such that S > 0, since 7 is to be 
negative. 

Note that since S is always positive for S, l z < 1 and 
for H = 0, the DL solutions still exist in the absence 
of H . This implies that quantum effects will be rele- 
vant in dense plasma environments where electrons and 
positrons are degenerate (e.g., in magnetars as well in the 
next generation laser solid-density plasma interaction ex- 
periments) for which classical fluid model fails to describe 
the plasma dynamics. Basically, the quantum parameter 
H (> H c ) restricts here the particle density to be of 
the order of 10 34 m -3 or lower (since higher values of H 
corresponds to lower density regimes) in order that the 
ion-acoustic DLs exist. Furthermore, the DLs are com- 
pressive or rarefactive according as 5 *s 0.001. Inspecting 
Eq. (|9|) one finds that while the width decreases, the 
amplitude of the stationary DLs increases with a slight 
increase of the density ratio 8. Also, the width increases, 
with the frequency lo c and the obliqueness l z . 

Next, we numerically investigate Eq. © using 
Runge-Kutta scheme with an initial condition </> (£, 0) = 

- (00/2) tanh ^-7/8 {(j>o/N) £ and with 

N = 3, 8 ~ 0.001, H - 0.22 
the coefficients in Eq. are A = —25.5, /1 = 0.17 and 
7 = —2.6 x 10 3 . Note here that though the initial condi- 
tion looks very similar to the asymptotic solution ([5]), but 
the constant N appears in different manner (For detailed 
discussion see, e.g., Ref. [33j]). We use 1000 grid points 
with the system scale size = 100, and choose the pulse 
size L p (~ 1.2) to be less than in order that the shock 
solutions exist. The shock profiles are shown in Fig. 1 
giving a spatial scaling constant (3.7 Xf) 1 for the ini- 
tial profile and (1.2Af) _1 for the asymptotic shock. The 
solitary waves in front of the shocks are ordered so that 
one with the maximum amplitude is nearest the shock. 
Initially the shock profile steepens, as time goes on it 
reaches a steady state with two solitary wave fronts (in 
the case of N = 3). 

In conclusion, the formation of forward propagating 
DL shocks is possible in a strongly magnetized quantum 
e-p-i plasma. Such shocks composed of the compressive 
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FIG. 1: (Color online) The profiles of the DL shocks for three 
different times. The initial profile (r = 0) steepens (dashed 
line), and simultaneously two solitary waves form in front of 
the shock (see the solid line for r = 150 and the dotted line 
for r = 200). 



as well as rarefactive slow- wave fronts propagate with the 
speed less than that of the ion-acoustic speed. Further- 
more, the DLs exist in dense plasma environments with 
particle density of the order of 10 34 m~ 3 or lower when 
the background electron population is much larger than 
that of positrons, and electron-positron annihilation is 
negligible [30]. Existence of such DL shocks may have 
significant role for the particle acceleration in the next 
generation laser solid-density plasma interaction experi- 
ments [j| as well as in compact astrophysical objects as 
evident from the recent observations in magnetized white 
dwarfs [l4j]- However, conclusive evidence needs further 
investigation in this area. In this way one may extend 
our investigation by considering the relativistic as well as 
the spin quantum effects [TiTj in a quantum MHD model, 
which, we hope, will give better understanding for the 
existence and properties of such novel DL structures. Fu- 
ture research is also expected to reveal other interesting 
applications of the shock solutions demonstrated here. 
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